1 



PSI-PR-02-09 



Electromagnetic radiative corrections to pionic beta decay rr 

H. Pichl 

Paul Scherrer Institut, Theory Group, CH-5232 Villigen PSI, Switzerland 



+ 



7r°e+z/ e * 



(N 
O 
O 
(N 

Oh' 

<D 

in 

(N 



> 
m 

On 
(N 
<3\ 
O 
(N 
O 

a: 

Oh! 



X 



Pionic beta decay 7r + — » 7r°e + ^ e is analyzed in chiral perturbation theory with virtual photons and leptons. 
All electromagnetic corrections up to 0(e 2 p 2 ) are taken into account. Theoretical results are confronted with 
preliminary data from a PSI measurement and a value for the CKM matrix element \V u d\ is given. Although 
the precision is presently still below the one of existing determinations of \V u d\, an analysis of pionic beta decay, 
based on a systematic treatment within a low-energy effective field theory, may become a useful alternative. 



1. Introduction 



One of the cornerstones of the Standard Model 
is the CKM matrix V. Its matrix elements V,j 
determine the strengths of transitions between 
quark flavours i and j. The Vy are therefore of 
eminent importance, and their precise determina- 
tion is a major task of today's particle physics. 

Unitarity implies a series of relations among the 
matrix elements. In particular, one expects the 
first-row CKM matrix elements to satisfy 
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|V M& | 2 = 1, 



(1) 



where mixing among the first-generation quarks 
u and d is strongest, and mixing between the first 
and third generation is very small and can be ne- 
glected at the present precision |Q. 

The strength of an s to u transition is given by 
V us which is measured in decays. Taking into 
account a recent analysis of radiative corrections 
in Km decays in chiral perturbation theory 
the Particle Data Group (PDG) leaves the central 
value of | Vis | unchanged but increases the error 
and gives finally |V US | = 0.2196 ± 0.0026 |J. 

The most precise up-to-date knowledge of \ V u d\ 
comes either from analyses of nuclear beta decays, 
in particular super-allowed Fermi decays, or (not 
as precise) from neutron decay. 

From beta decays of nuclei one finds |V«;| = 
0.9740 ± 0.0005, where the uncertainty depends 
on structure-dependent radiative corrections JD. 



*Talk presented at QCD 02, Montpellier, 2-9 July 2002. 
This article is based on work done together with V. 
Cirigliano, M. Knecht and H. Ncufcld. 



Interactions with the nuclear medium could give 
the quarks different effective masses, which might 
lead to an enhanced value of \V uc i\- The PDG 
compensates for this possibility by doubling the 
aforementioned error |Q. 

An extraction of | V u d | from neutron decay ben- 
efits on one hand from fewer theoretical problems, 
but suffers at the same time from the fact that 
not only the neutron lifetime but also the ratio 
of axial- vector and vector couplings g& / gv comes 
into play. Averaging over neutron decay experi- 
ments (see [Q and references therein), the PDG 
gives \V u d\ = 0.9725 ± 0.0013, and a final com- 
bined value of \V ud \ = 0.9734 ± 0.0008 with a 
remarkably small error is found 0. 

With these values one finds a 2.2er deviation 
from unitarity in (@): | V ud \ 2 + \V US | 2 + 1 V ub \ 2 - 1 = 
-0.0042 ± 0.0019. Adopting a conservative line 
of reasoning, one may speculate if the given error 
of | V u d | might be too optimistic, and look for a 
cleaner theoretical alternative to determine |Vid|- 

Pionic beta (tt/3) decay ir + — > Tr°e + v e combines 
the advantages of the traditional approaches, i.e., 
it is a pure vector transition without any compli- 
cations due to nuclear structure. However, practi- 
cal utilization of nf3 decay is severely complicated 
by its branching ratio of BR^ ~ 1 x 10~ 8 [||J]. 

In spite of the rareness of the decay, a measure- 
ment of BR^ is being performed at PSI, aiming 
to achieve a precision of 0.5% in a first stage [Q. 
Of course, this first accuracy will not be enough 
to compete with nuclear or neutron decays, but 
a third independent determination of \V u d\ may 
certainly be useful. 
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At this level of accuracy radiative corrections 
have to be included in the theoretical analysis, 
and we employ chiral perturbation theory as the 
ideal framework for low-energy hadron physics^]. 
This article is based on (5) which in turn is based 
on recent work on K13 decays Q . 

2. Kinematics and decay amplitude 

Neglecting radiative corrections for the mo- 
ment, the decay amplitude of ir/3 decay 



■X + (P+) -> 7T°(Po) e + (p e ) V e {Pu) 

is written in terms of two form factors : 



(2) 



m = g f v;j, 



(0) 



where = u(p v )^^(l — 75)w(p e ) is the lepton 
current. Gf is Fermi's coupling constant. 

At leading order, /| 0) = 1 and f W = 0. At 
next-to-leading order, the form factors develop a 
dependence on the variable t = (p+ — Po) 2 - 

/_ turns out to be proportional to the dif- 
ference of squared pion masses (m + 

Additionally, /i "* terms in the squared decay 
amplitude come always along with a factor of 
r e = m 2 ./m 2 . + ~ 1.35 x 10~ 5 . These observations 

justify to neglect /_ in the amplitude. 
The decay rate T^p for 7r/3 decay, 



= / dydz p( Q) (y,z) 
Jv 



(4) 



is expressed in terms of a phase space density p^ 
which depends on two variables y and z defined 
in the c.m.s. by 

2E„ . 2E„a 

V 



and 



777,,+ 



(5) 



T> denotes the kinematically allowed area of (y, z) 
The phase space density itself is given by 

pM(y,z)=Afx\fW(t)\ 2 AW(y,z), 



(6) 
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where A^°\y : z) is a kinematical density 

A {0 \y, z) = [4 (-1 + z + 2y -yz - y 2 ) - 
+ r e (Ay + 3z - 3 + r - r e )] , 

1 Jaus also Q addressed the problem of rad. corrections 
to 7T/3 decay in a constituent quark model. 



(7) 



and M in @ and ro in (Q) are defined by 

., \Vu d \ 2 G F ml + 

64tt 3 1 r =m w0 /m„+. (8) 

At 0(p ), meson loops and counterterms from 
the chiral Lagrangian £ 4 generate the afore- 
mentioned t-dependencef| (t) = 1 + Sf^ (t) , 
1 



F 2 



[2tL r Q (p) + 2h r ^ (t,p) 



+h r K+K0 (t,p)}. (9) 

The functions h r P Q (t, p) comprise one- loop cor- 
rections from the particles P and Q |?J . Both the 
renormalized low-energy coupling (LEC) Lg(p) 
and h P Q(t, p) depend individually on the ren. 

scale, but Sfi. (t) is independent of p. 

3. Radiative corrections 

Radiative corrections manifest themselves in 
virtual photon exchange, new electromagnetic 
counterterms, and in the emission of real photons. 

Photon exchange generates one-loop diagrams 
as a result of which a new dynamical variable u = 
(p+ —pe) 2 and infrared divergences appear in the 
amplitude M. Counterterms from Lagrangians 



£„2„2 and C\, 



contribute three new LECs 



X u X£(p), and K{ 2 (p). 

Virtual radiative corrections induce a change 
of the form factor from f+\t) to f + (t,u,X) = 

f+ (t) + <$/+ m ( u >-M> where A plays the role of a 
fictitious photon mass to regularize infrared di- 
vergences. 

Only showing local contributions explicitly and 
omitting remaining one-loop contributions, the 
radiative corrections of 0(a) read 



Sf. 



2 1 
i..,al = 47ra<j -~Xi - -Xl(p) 



2K{ 2 {p) 



32tt 2 



3 + In ■ 



3 In- 



(10) 



As before, <5/?. m (t4, A) is independent of the renor- 
malization scalc^J. 

2 O(e 2 p ) corrections are absorbed in the meson masses. 
3 The omitted one-loop part Sf""' (u, A) containing in par- 
ticular the infrared divergences can be found in |5l . 
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As it was shown in §|| , it is useful to rearrange 
the complete form factor in the following way: 

/+(*, u, A) = F+(t) x [1 + 5ffi(u, A)] , (11) 

by factoring-out the ^-dependence in F+(t) := 
/+°^ (t) + Sf+i ocai and isolating both infrared di- 
vergences and M-dependence. This allows us to 
write the rate in the presence of virtual radia- 
tive corrections in combination with real photon 
emission formally in exactly the same way as in 
(Q), (6), i.e., in terms of a generalized form factor 
F+(t) and a generalized infrared finite kinemati- 
cal density. 

The treatment of real photon emission and can- 
cellation of infrared divergences was discussed in 
detail in Working at order a, it is suffi- 

cient to consider the emission of only one real 
photon. We adopt a scheme proposed by Gins- 
berg Q and integrate over the entire phase space 
of undetected particles, i.e., the photon and the 
neutrino, but restrict the (y, z) space to the non- 
radiative three-particle decay case. Furthermore, 
a new variable x = (p„ + p 7 ) 2 is defined and in- 
tegrated over. 

The phase space density for the radiative decay 
7r + — > 7r°e + ^ e 7 is given by 



2 13 7T 6 



fix 



\ 2 



E v E 1 



+ - Po - Pe - P„ - P 7 ) ^ l-M Y , (12) 

pol 

where M 1 denotes the leading-order amplitude 
of the radiative decay, respectively. Upon com- 
bining the phase space densities of the mother 
process and the associated radiative decay, one 
arrives at an infrared finite phase space density 

p(y,z) = pW(y,z) + pr(y i z) (13) 

which we may write, as already anticipated, as 



p{y,z)=Afx \F+(t)\ 2 xA(y,z). 



(14) 



in terms of a generalized kinematical density 
A(y,z) and the form factor F + (t). Apart from 
electromagnetic local contributions, radiative cor- 
rections from virtual photon exchange ( |l0| ) and 
real photon emission ([l2|) are entirely included in 
A{y,z) ' 



4. Results and conclusions 

Before presenting our results, a few comments 
seem advisable. In general, low-energy couplings 
of chiral Lagrangians parametrize short-distance 
physics that is not explicitly dealt with in the 
effective theory. 

According to |0,0|, it turns out that all 
semileptonic charged current amplitudes are af- 
fected by universal short-distance physics correc- 
tions when expressed in terms of the muon de- 
cay constant. Relating therefore Gf with the 
muon lifetime, it is possible to trace back some 
of the high-energy origin of the LEC Xq to these 
universal corrections Hence, we may write 
= A| D + Jq(fj), and define e 2 A 6 SD = 
1 — S QW (m p ,mz) to contain the short-distance 
enhancement factor S ev ,(rn p ,mz) = 1.0232 from 
The remaining part Xg(/i) is now expected 
to be of the order of ~ 10 -3 . 

Extracting the short-distance physics from 
p(y,z) and expanding until 0(t), we may write 
the infrared safe inclusive decay rate as 

r ffj 8( 7 ) =MS cw {m pi m z ) x \F+{0)\ 2 I(\+). (15) 

The linear approximation works extremely well, 
since m\<t < (m^+ - m^o) 2 ~ 21.1 MeV 2 . 



= / dydzA(y,z) 



in- 



(16) 



contains the slope parameter A + from the expan- 
sion in t. 

To estimate F + (0), we use the 'classical' value 
of Ll(m p ) = (6.9 ±0.7) x 10~ 3 , quoted in §, and 
take Moussallam's estimate |l3| of i^[ 2 ( TO p) = 
(-4.0 ±0.5) x 10" 3 . 

Since essentially nothing^ is known about the 
Xi, we resort to naive dimensional analysis and 
use |Ai|, \Xl{m p )\ < 6.3 x 10~ 3 . 

Finally, F+(0) is found to be 1.0046 ± 0.0005. 
The error is completely due to uncertainties in 
the electromagnetic LECs and is extremely small. 
The slope parameter's uncertainty comes from 
the error of Lg(m p ), and we obtain A+ = (0.037± 
0.003). 



4 In a combined analysis of K^t an d K® 3 decays it is pos- 
sible to extract a value for Xi |14j. 
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The phase space factor is calculated to 

be J(A+) ~ 7.3832 x 10~ 8 with a tiny error. 

We arrive at the following results. Radiative 
corrections turn out to be at the < 1% level. The 
corrected phase space /(A+) enhances the rate by 
- 0.1%, the form factor F+(0) by ~ 0.9%. 

Using the PDG's value of \V ud \ = 0.9734 ± 
0.0008, we arrive at a 7r/3 branching ratio of 



BR n p = (1.0376 ±0.003) x 10" 



(17) 



Turning the procedure upside down, we use the 
7T/3 branching ratio together with our calculation 
of radiative corrections to extract 



\V ud \ = 9600.8 x 



l*+(0)| ' 

The error assigned to our value of \V u d\ is 



A|Kd| = |Kd| x 



± AF+(0) 1 ABR^ (7) 



F+(0) 



2 BR 



7T/3( 7 ) 



(18) 



(19) 



For the PDG's value of the branching ratio 
BRpdg = (1-025 ± 0.034) x 10~ 8 , we get 



w, 



Md th 1 



0.9675 ± 0.0005(th) ± 0.017(cxp). (20) 



Obviously, the theory is very well under control 
and it is the experiment that limits the precision. 

The situation is much improved if one uses 
new PSI data on ir/3 decay With the re- 

cent, however still preliminary, branching ratio 
of BRpgi = (1.044 ± 0.016) x 10~ 8 , we find 



I 



th2 



0.9765 ± 0.0005(th) ± 0.0075(exp).(21) 



From the theorist's point of view, we wish to 
emphasize that the inclusion of radiative correc- 
tions is extremely straightforward and particu- 
larly simple. The remarkably small theoretical 
error can practically entirely be attributed to the 
low-energy couplings X\ and Xg(m p ). 

Moreover, the approach is completely equiva- 
lent to the one employed in decays ]2j , which 
puts both processes on an equally sound footing. 
Note finally that the central value of \V u d\ in ( ]2l| ) 
matches the unitarity constraint (|l|) very well, but 
the error is, of course, much too large to allow for 
more stringent final conclusions. 

It is fair to say that the present experimental 
precision of 7r/3 decay cannot yet compete with 



nuclear beta decays, and even if the final pre- 
cision will drastically increase it remains rather 
unlikely to reach a competitive accuracy. Nev- 
ertheless, 7T/3 decay could become a valuable al- 
ternative source for the determination of \V ud \, in 
particular if previous error estimates turned out 
to have been too optimistic. 
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